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Abstract. We perform a global analysis of transverse momentum distributions in Drell-Yan pair 
and Z boson production in order to investigate universality of nonperturbative contributions to the 
Collins-Soper-Sterman resummed form factor. Our fit made in an improved nonperturbative model 
suggests that the nonperturbative contributions follow universal nearly-linear dependence on the 
logarithm of the heavy boson invariant mass Q, which closely agrees with an estimate from the 
infrared renormalon analysis. 

ON ' 

^ . Transverse momentum distributions of heavy Drell-Yan lepton pairs, W, or Z bosons 

produced in hadron-hadron collisions present an interesting example of factorization 
for multi-scale observables. If the transverse momentum qj of the electroweak boson 
is much smaller than its invariant mass Q, da jdqj at an n-th order of perturbation 
theory includes large contributions of the type a"ln m (qj/Q 2 )/qj< (m = 0, 1 . . .In — 1), 
which must be summed through all orders of Ot s to reliably predict the cross section [1]. 
Such resummation is realized in the Collins-Soper-Sterman (CSS) formalism [2], which 
describes soft and collinear QCD radiation in a wide range of energies by introducing a 
resummed form factor W(b) in impact parameter (b) space. 

While the short-distance contributions (b < 1 GeV -1 ) to the CSS form factor W(b) 
can be calculated in perturbative QCD, long-distance nonperturbative contributions from 
b > 1 GeV 1 are not yet fully computable, even though their basic form can be deduced 
from the infrared renormalon analysis [3]. The factorization theorem behind the CSS 
formalism predicts that the nonperturbative contributions are universal in unpolarized 
Drell-Yan-like and semi-inclusive DIS processes. Consequently the function ^Mp{b, Q) 
that describes the nonperturbative terms can be constrained in a global fit to the hadronic 
qr data, just as the fcj -integrated parton densities are constrained with the help of 
inclusive scattering data. ^Np{b,Q) must be known precisely in order to successfully 
measure the W boson mass, because uncertainties in J?Np(b, Q) may affect the measured 
value of Myy at the level comparable to the targeted accuracy of the measurement, 
8M\y « 30 MeV at the Tevatron and 15 MeV at the LHC. It is therefore interesting to 
investigate if ^Np{b, Q) found in the qj fit is consistent with the universality hypothesis, 
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and whether its preferred form is compatible with the renormalon analysis. 

These issues were explored recently in Ref. [4], where a global analysis of q T data 
from fixed-target Drell-Yan pair production and Tevatron Z boson production was per- 
formed in the context of an improved model for the nonperturbative contributions. Al- 
though ^^p{b, Q) primarily parametrizes the "power-suppressed" terms, i.e., terms pro- 
portional to positive powers of b, its form found in the fit is correlated with the assumed 
behavior of the leading-power terms (logarithmic in b terms) at b < 2 GeV 1 . The exact 
behavior of W(b) at b > 2 GeV -1 is of reduced importance, as W(b) is strongly sup- 
pressed at such b. For these reasons, we closely followed the procedure of the previous 
global qr analysis [5], while paying close attention to the model of the leading-power 
terms at perturbative and moderately nonperturbative transverse distances, b < 2 GeV 1 . 

The large-Z? contributions were introduced by using the b* model [2], as 

W(b)=W pert (b,)e-^ b ' Q l (1) 

Here W perl (Z?*) is the perturbative part of W (b) , i. e. , its leading -power part evaluated at a 

finite order of a s . W pert (b*) depends on the variable b^ = b/(l + b 2 /b 2 nax ) 1 / 2 and serves 
as an approximation for all leading-power terms. Its shape is varied at all b by adjusting 
a single parameter b max . The b* model with a relatively low b max = 0.5 GeV -1 was 
a choice of the previous qj fits [5, 6]. However, it is natural to consider b max above 
1 GeV -1 in order to avoid ad hoc modifications of W pert (b) in the b region where 
perturbation theory is still applicable. In Ref. [4], we proposed a modification in the 
b* model that allowed us to increase b max at least up to ps 3 GeV -1 , while preserving 
correct resummation of the large logarithms at small b and numerical stability of the 
Fourier-Bessel transform. If a very large b max comparable to I/Aqcd is taken, Wip(b) 
essentially coincides with W pert (b), extrapolated to large b by using the known, although 
not always reliable, dependence of W pert (b) on \nb. Hence, the new prescription can be 
also used to test viability of extrapolation of W pert (b) to large b, reminiscent of similar 
extrapolations introduced in the alternative models [7, 8]. 

Following the renormalon analysis and Ref. [5], we assumed a Gaussian form of the 
nonperturbative function, <^Np{b,Q) = a(Q)b 2 , with 

a(Q) = fli + a 2 ln [2/(3.2 GeV)] + a 3 ln [lOCki^] . (2) 

The dependence of J^jyp on In Q is a consequence of renormalization-group invariance 
of the soft-gluon radiation. The coefficient of the In Q term has been related to 
the vacuum average of the Wilson loop operator and evaluated within lattice QCD as 
0- 19lo'o9 GeV 2 [9]. To see if the universal Gaussian behavior is consistent with the data, 
we first examined the values of a (Q) that are independently preferred by each bin of Q in 
5 examined experimental data sets. Fig. 1(a) shows the best-fit values of a(Q) obtained 
in independent fits to the data in each bin of Q for b max = 1.5 GeV -1 . The best-fit a{Q) 
follow a nearly linear dependence on In Q, and the slope at = da(Q) /d (In Q) is close to 
the renormalon analysis expectation of 0.19 GeV 2 [9]. Such nearly linear behavior of 
a (Q) is observed in the entire range b max =1—2 GeV -1 , and it less pronounced at b max 
outside of the interval 1-2 GeV -1 . Since the best-fit a(Q) in each Q bin are essentially 
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Figure 1. (a) The best-fit values of a(Q) obtained in independent scans of x 2 f° r the contributing 
experiments. The vertical error bars correspond to the increase of % 2 by unity above its minimum in 
each Q bin. The slope of the line is equal to the central-value prediction from the renormalon analysis [9]. 
(b) The best-fit % 2 and coefficients a\, a^, and a$ in ^Np{b,Q) for different values of b max . The size of 
the symbols approximately corresponds to lc errors for the shown parameters. 



independent, we conclude that the data support the universality of ^mp, when b max lies 
in the range 1—2 GeV -1 . In addition, each experimental data set individually prefers a 
nearly quadratic dependence on b, J^Vp = a(Q)b 2 ^ , with |j8| < 0.5 in all experiments. 

Next, we performed a simultaneous fit of our model to all the data. Fig. 1(b) shows 
the dependence of the best-fit x 2 a\, a2, and a-$ on b max . As b max is increased above 
0.5 GeV -1 assumed in the studies [5, 6], X 2 rapidly decreases, becomes relatively flat at 
bmax =1—2 GeV -1 , and grows again at b max > 2 GeV -1 . The global minimum of x 2 
is reached at b max ^1.5 GeV -1 , where all data sets are described equally well, without 
major tensions among the five experiments. The magnitudes of a\, ci2, and a?, are reduced 
when bmax increases from 0.5 to 1.5 GeV -1 . In the whole range 1 < b max < 2 GeV -1 , 
«2 agrees with the renormalon analysis estimate. The coefficient aj, which parametrizes 
deviations from the linear In Q dependence, is considerably smaller (< 0.05) than both 
a\ and ci2 (~ 0.2). This behavior supports the conjecture in [7] that as is small if the 
exact form of W pert (b) is maximally preserved. 

The preference for the values of b ma x between 1 and 2 GeV -1 indicates, first, that 
the data do favor the extension of the b range where all leading-power terms are ap- 



proximated by their finite-order expression W pert (b). In Z boson production, this region 
extends up to 3 — 4 GeV 1 as a consequence of the strong suppression of the large-£ tail 
by the Sudakov exponent. The fit to the Z data is actually independent of b max within 
the experimental uncertainties for b max > 1 GeV -1 . In the low-<2 Drell-Yan process, the 
continuation of bW pert (b) far beyond b pa 1 GeV -1 is disfavored because of large higher- 
order corrections to bW pert (b) at b around 1.5 GeV -1 . To summarize, the extrapolation 
of W pert (b) to b > 1.5 GeV -1 is disfavored by the low-<2 data sets, if a purely Gaussian 
form of J^Xrp is assumed. The Gaussian approximation is adequate, on the other hand, 
in the b* model with b max in the range 1 — 2 GeV 1 . 

In Z boson production, our best-fit a(M z ) = 0.85 ±0.10 GeV 2 agrees with 0.8 GeV 2 
found in the extrapolation-based models [7, 8], and it is about a third of 2.7 GeV 
predicted by the BLNY parametrization. In the low-<2 Drell-Yan case, our a(Q) = 
0.2 - 0.4 GeV 2 is close to the average (a) = 0.19 - 0.28 GeV 2 in four Q bins of the 
E288 and E605 data found in the model [7]. To describe the low-<2 data, Ref. [7] allowed 
a large discontinuity in the first derivative of W(b) at b equal to the separation parameter 

bmax = 0.3 — 0.5 GeV -1 , where switching from the exact W pert (b) to its extrapolated 
form occurs. In the revised b* model, such discontinuity does not happen, and WLp{b) is 
closer to the exact W perl (b) in a wider b range than in Ref. [7]. 

The best- fit parameters in found in the new model are quoted in Ref. [4]. The 
global fit places stricter constraints on J^vp at Q = Mz than the Tevatron Run-1 Z 
data alone. Theoretical uncertainties from a variety of sources may be substantial in 
the low-<2 Drell-Yan process, which is indicated, in particular, by the dependence of 
the agreement with the low-<2 data on an arbitrary factorization scale C3 in W pert (b). 
The low-<2 uncertainties do not substantially affect predictions at the electroweak scale. 
The &{c^) corrections and scale dependence are smaller in W and Z production, and, 
in addition, the term ci2lnQ, which arises from the soft factor y(b,Q) and dominates 
<^np at Q = M z , shows little variation with C3. Consequently, the revised b* model with 
bmax ~ 1-5 GeV -1 increases our confidence in the transverse momentum resummation 
at electroweak scales by exposing the soft-gluon origin and universality of the dominant 
nonperturbative contributions at collider energies. 
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